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AN  EXAMINATION  OF  MODELS  OP  RELAXATION  IN  COMPLEX  SYSTEMS 
L  Continuous  Time  Random  Walk  (CTRW)  Models 

A.  Introduction 

We  Initiate  here  a  detailed  examination  of  existing  models  of  relaxation 
in  glass-forming  materials.  In  general^  our  considerations  will  include:  (1) 
the  robustness  of  the  models;  (2)  technical  procedures  employed;  (3)  physical 
tests  of  asstmptions;  (4)  predictability.  A  model  is  robust  if  it  is  insensi¬ 
tive  to  untestable  or  omitted  details  while  being  consistent  with  tests  of 
premises  and  established  facts.  The  predictability  of  a  model  is  determined 
by  its  capability  to  generate  new  predictions  that  are  actually  verified 
empirically. 

In  this  report  we  consider  relaxation  models  [1,2]  based  on  the  continuous 
time  random  walk  (CTRW)  formalism  [3].  In  these  models,  mobile  random  walkers 
identified  as  defects  move  on  a  lattice  in  a  sequence  of  steps  and  pauses,  and 
cause  relaxation  of  the  entities  of  interest  which  are  located  at  lattice 
points.  The  entities  of  interest  will  be  designated  the  relaxors.  The  basic 
function  entering  into  these  models  is  the  pausing  time  distribution  <|;(t) 
which  must  be  specified  either  by  assumption  or  by  a  further  model  within  a 
model.  Although  the  nature  of  the  defects  and  the  physical  mechanism  of 
defect-relaxor  interaction  is  left  unspecified,  the  models  based  on  CTRW  do 
include  an  assumption  about  whether  N,  the  number  of  defects  per  relaxor,  is 
large  or  not. 

As  will  be  discussed  below,  for  a  given  (ii(t),  different  relaxation  func¬ 
tions  are  obtained  depending  on  which  of  the  above  assumptions  about  N  are 
made  [4,5,6].  Since  typically  no  information  is  provided  on  the  identities  of 
Manuw-ipt  approved  October  31,  1985. 


the  defect  and  the  relaxor,  either  assumption  oust  be  considered  ad  hoc  physl-* 
cally.  However,  there  are  physical  systems  which  are  also  glass-formers  and 
have  mobile  entitles  which  can  be  considered  as  random  walkers  In  a  CTRW 
framework.  Such  physical  situations  provide  a  useful  testing  ground  for  the 
validity  of  the  assumptions  used  In  relaxation  models  based  on  CTRW.  Then  the 
CTRW  models  can  be  evaluated  according  to  the  considerations  mentioned  above: 
(1)  Their  robustness  can  be  tested  on  the  basis  of  the  sensitivity  of  the 
dependence  of  the  calculated  relaxation  function  on  the  asstnptlon  made  for 
form  of  the  pausing  time  distribution  (l>(t),  and  also  on  the  assumption  made 
for  N,  the  ntimber  of  defects  per  relaxor.  (2)  Their  technical  soundness  can 
be  assessed  by  careful  evaluation  of  the  assumptions  used  and  modelling  pro¬ 
posed  for  (l;(t).  (3)  The  general  assumptions  can  be  tested  physically  by  their 
consistency  with  observed  phenomena  In  physical  situations  that  Involve  random 
walkers.  (4)  Predictability  can  be  tested  by  either  the  comparison  of  pre¬ 
dicted  phenomena  with  observed  results  or  the  absence  of  descriptive  power  In 
dealing  with  established  phenomena. 

The  relaxation  model  based  on  CTRW  proposed  by  Shleslnger  and  Montroll 
(SM)  [1]  which  Is  a  generalization  of  earlier  defect  diffusion  models  of 
Glarum  and  others  [7]  Is  representative  of  the  class  of  CTRW  relaxation 
models.  We  will  focus  attention  on  the  SM  model.  A  representative  example  of 
modelling  of  (|;(t)  has  been  provided  by  Bendler  and  Shleslnger  (BS)  [2]  ,  and 
this  work  will  also  be  reviewed  In  detail. 

B .  Robustness 

As  mentioned  above,  a  model  of  dielectric  relaxation  based  on  CTRW  has 
been  promulgated  by  Shleslnger  and  Montroll  (SM)  [1] .  This  model  is  a  gener¬ 
alization  of  the  defect  diffusion  models  of  Glarum  and  others  [7]  which  assume 
that  the  physical  system  contains  two  types  of  physical  species:  frozen 


dipoles  and  mobile  defects.  On  reaching  the  site  of  a  frozen  dipole,  a  mobile 
defect  causes  the  dipole  to  relax.  That  Is,  a  contact  Interaction  between 
dipole  and  dipole  Is  assumed  but  no  physical  mechanism  Is  specified. 

In  the  SM  model,  a  frozen  dipole  Is  at  the  origin  of  a  lattice  and,  at 
t*0,  the  defects  lie  with  equal  probability  at  any  lattice  point.  Diffusion 
of  the  defects  Is  executed  as  a  continuous-time  random  walk  [3]  composed  of  a 
sequence  of  steps  and  pauses.  The  pausing  time  distribution  function  (];(t)  is 
assumed  to  have  a  power  law  form  at  sufficiently  long  time 

(^(t)~At-l“a,  0<a<l,  t->.«  (1) 

The  complete  distribution  function  must  be  normalizable 

^/-dtcKt)  .  1.  (2) 

Obviously,  the  form  of  Eq.  (1)  cannot  be  extrapolated  back  to  short  times  and 
still  satisfy  the  normalization  requirement  of  Eq.  (2).  The  form  of  (i)(t)  at 
short  times  has  not  been  specified  In  the  model.  As  will  be  detailed  below, 
the  model  [1]  then  commits  to  the  assumption  that  N,  the  number  of  defects  per 
dipole.  Is  large.  The  defects  and  the  dipole  are  not  assigned  explicit  physi¬ 
cal  identities.  Therefore,  any  assumption  concerning  the  relative  number  of 
defects  (e.g.  large,  equal,  etc.)  to  dipoles  or  relaxors  does  not  have  a 
direct  physical  basis.  Sensitivity  to  the  assumption  about  relative  number  of 
defects  that  will  be  shown  to  exist  reveals  one  lack  of  robustness  In  CTRW 
models. 

SM  follow  Tachiya  [6]  in  the  derivation  of  the  equation  for  the  survival 
probability  of  dipoles.  Let  Njj  be  the  number  of  frozen  dipoles  in  a  polymeric 
system  or  a  glass  of  volume  Q.  The  system  is  divided  into  Njj  identical  unit 
cells  each  of  volume  Q^^Q/Nj,  and  each  containing  one  frozen  dipole  and  N 
defects.  Each  dipole  is  fixed  at  the  origin  of  Its  unit  cell  and,  excluding 
the  origin,  there  are  V  lattice  points  for  each  unit  cell.  The  N  defects  in 


each  ualt  cell  are  treated  as  Independent  and  are  totally  randomly  distributed 

on  the  V  lattice  sites.  Periodic  boundary  conditions  proposed  by  Shleslnger 

[5]  are  used  In  each  unit  cell  such  that  If  a  defect  leaves  the  cell.  It  is 

effectively  replaced  so  that  the  proper  defect  number  N  In  each  unit  cell  Is 

maintained.  Let  f(l,t)  be  the  probability  density  that  a  defect  originally  at 

a  lattice  point  1  reaches  the  origin  for  the  first  time  at  time  t.  Then 

[l-^/^dt'f(l,t' ) ]  Is  the  probability  that  this  defect  has  not  reached  the 

origin  In  the  time  Interval  (o,t).  The  defects  are  treated  as  Independent  so 

so  there  is  a  tacit  assumption  that  they  are  non-interacting  and  there  are  no 

correlations  among  the  defects.  The  quantity  n  fl-  /  dt' f( 1. ,t' ) 1  Is  the 

1=1  0  1 

probability  that  none  of  the  N  defects  has  reached  the  origin  In  the  time 
interval  (o,t).  Averaging  over  the  initial  uniform  distribution  of  defect 
positions  gives  an  expression  for  the  survival  probability  of  frozen  dipoles 
(the  relaxation  function  (|)(t) )  as 

<{,(t)  -  [1-  Z  /dt’f(Jl,t’)f-[l-  4  At'I(t*)f 


(3) 


Here  I(t)=  Z  f(l,t)  is  the  flux  of  defects  into  the  dipole.  This  result  has 
been  derived  as  Equation  (20)  in  Reference  [5]  and  equation  (3.5)  in  Reference 
[6] ,  and  written  as  Eq.  (8)  in  the  work  by  SM.  Now  we  come  to  a  step  in  which 
a  critical  assumption  about  the  size  of  N  has  to  be  made. 

If  the  number  N  is  large,  the  right  hand  side  of  Eq.  (3)  can  be  approxi¬ 
mated  by  an  exponential  function.  With  constant  defect  concentration  c=N/V  as 
N-^a>,  V->oo,  the  dipole  relaxation  becomes 

4)(t)  «  exp[-c  /^dt'I(t')] 


SM  complete  their  dielectric  decay  model  by  relating  the  flux  I(t)  to  the 
pausing  time  distribution  clj(t)  of  the  defects.  The  relation  between  the 
Laplace  transforms  of  I(t)  and  (j;(t),  namely  I(u)  and  ci;*(u),  turns  out  to  be 


'  '  [1-(1;*(U)]P[0,(|,*(U)  j  »  (5) 
where  the  specific  form  of  the  function  P(o,z)  depends  on  the  type  and 

dimension  of  the  lattice.  SM  have  only  evaluated  I(u)  for  small  u  which 
corresponds  to  the  long  time  regime.  They  cannot  address  the  short  time 
regime  because  only  the  long  time  behavior  of  4;(t)  has  been  specified.  For  a 
simple  cubic  lattice  In  three  dimensions,  the  SM  method  of  evaluating  I(u) 
leads  to  l(u)“0. 659u“o/A  for  small  u.  This  Is  given  by  Eq.  (27)  In  [1].  This 


leads  to 


I(t)~t®'”^  ,  t->® 


(5SM) 


However,  In  an  earlier  work  on  the  electron  scavenging  problem,  Shleslnger  [5] 
used  an  expression  which  Is  completely  equivalent  to  Eq.  (3)  and  where 
(|j(t)~t"l“a  as  t-^ca.  In  this  earlier  work,  the  expression  for  I(u)  for  small  u 
was  evaluated  In  a  different  way  and  Shleslnger  obtained  I(u)“l-Bua  for  small 
u.  This  Is  given  by  Eq.  (17c)  In  [5]  where  F(u)  In  [5]  Is  proportional  to 
I(u)  here.  Shleslnger  uses  this  to  obtain  the  relation  11m  I(t)~llm  (|;(t),  so 

t-*'08  t-V® 

that 

I(t)~t"a"l  ,  t-*-®  (5S) 


The  result  for  I(t)  as  t-*-®  Is  quite  different  depending  on  whether  I(u)  Is 
evaluated  at  small  u  by  the  method  of  SM  or  by  the  earlier  method  of  Shlesln¬ 
ger.  Since  both  these  asymptotic  expressions  begin  from  the  same  equations. 
It  Is  Incumbent  upon  these  authors  [1,5]  to  clarify  which  method  for  obtaining 
I(u)  for  small  u  Is  appropriate. 

If  the  SM  expression  Eq.  (5SM),  I(t)~t®“l  as  t-*’®,  Is  used  In  Eq.  (4),  we 


(t((t)~exp(-cita)  ,  (6) 

where  ci  ®  0.659c/ (r(a)r(l-a)A].  The  SM  model  thus  leads  to  the  stretched 
exponential  for  ij)  but  one  that  Is  only  valid  as  t-*-®,  i.e.  at  the  tall  end  of 


\  ■* 


>  V  •  V"-  W*’  V."-  V*-  «*•  *>"'•  V  •  w'-  V* 


the  stretched  exponential.  Note  that  if  the  earlier  Shlesinger  expression 
Eq.  (5S),  I(t)'”t“®"^  as  t><»,  had  been  used  in  Eq.  (4),  we  would  have  obtained 
(()(t)~e:q>(-c2t"“)  ,  t^<»  (6S) 

Instead  of  the  stretched  exponential  Eq.  (6). 

The  stretched  exponential  Eq.  (16)  for  both  short  and  long  times  was  first 
proposed  by  Kohlrausch  [8]  in  1847.  This  is  now  often  referred  to  as  the 
Kohlrausch-Williams  and  Watts  (KWW)  function  [9].  Incidentally,  the  asymcotic 
large  t  behavior  of  (t)(t)  as  given  by  relation  (6)  was  actually  derived  earlier 
by  Hamlll  and  Funabashi  [4]  using  the  same  tall  for  4)(t)  as  in  Eq.  (1)  and  the 
same  conditions  including  N»I  for  the  electron  scavenging  problem  [10]  .  How¬ 
ever,  Tachiya  [6]  has  shown  that  the  Hamlll- Funabashi  model  involving  diffu¬ 
sion  of  either  the  trapped  electron  or  the  acceptor  molecule  is  not  relevant 
for  explanation  of  the  unusual  kinetics  of  trapped  electron  decay.  Tachiya 
maintained  that  a  tunneling  mechanism  proposed  earlier  and  which  Invokes 
direct  electron  transfer  from  trap  to  acceptor  [6,10]  is  the  correct  one.  In 
any  event,  such  a  mechanism  can  be  used  to  explain  the  experimental  data  [10] . 

It  is  Important  to  emphasize  that  the  limited  result  of  an  asymptotic  KWW 
function  obtained  by  the  SM  model  is  even  further  weakened  if  N  is  finite  and 
not  actually  taken  to  the  limit  for  fixed  c“N/V.  For  finite  N,  we  must 
return  to  Eq.  (3).  If  the  tail  in  Eq.  (1)  is  valid  for  times  greater 

than  some  time  t*,  then  we  can  use  the  previous  result  of  SM,  Eq.  (SSM),  that 
I(t)at®"^  for  t>t*  to  express  (t)(t)  in  Eq.  (3)  as 

♦  (t)  - 

Here  z  *  c  dt'I(t' )-c^t*®.  Recall  that,  according  to  the  SM  model,  the 
o  ^ 

2  a  - 

quantity  in  the  brackets  [l-  ^  t“]  represents  the  probability  that  the 

I 

defect  has  not  reached  the  dipole  at  time  t.  For  t>t**  =  [(N-z)/ci]  this 
survival  probability  would  become  negative  so  that  the  formalism  cannot  be 


applied.  Therefore,  when  N  Is  finite,  Eq.  (7)  is  restricted  in  its  applicabi¬ 
lity  by  both  a  lower  and  an  upper  time  limit;  t*<t<t**.  This  is  an  inherent 
problem  for  the  SM  model  with  finite  N.  For  sufficiently  large  but  fixed  N, 
it  is  possible  for  Eq.  (7)  to  approximate  a  KWW  function  over  some  time  regime 
between  the  limits  t*  and  t**.  Tachiya  has  pointed  out  [6]  that,  whatever  the 
value  of  N,  the  approximation  of  [l-  to  exp(-x)  becomes  increasingly  worse 
as  x->-<».  Therefore,  for  fixed  large  N,  the  time  interval  over  which  Eq.  (7) 
approximates  a  KUW  function  depends  both  on  the  value  of  N  and  the  values  of 
t*  and  t**.  And  when  N  is  not  sufficiently  large,  Eq.  (7)  will  not  approxi¬ 
mate  a  KWW  function  at  any  time,  even  between  the  limits  t*  and  t**. 

An  interesting  special  case  of  Eq.  (7)  occurs  for  N»l,  i.e.  when  there  is 
only  one  defect  per  relaxor.  Then  the  relaxation  function  is  given  by 

4)(t)  “  [l-z-cito]  (7a) 

Just  as  Eq.  (7)  this  result  follows  from  Eq.  (5SM)  (Eq.  (28)  of  SM) .  An  alter¬ 
native  procedure  for  evaluating  I(t)  was  used  earlier  by  Shlesinger  [5]  which 
leads  to  Eq.  (5S)  (Eq.  (22)  of  Ref.  [5]).  Shlesinger  rewrote  Eq.  (3)  as  ii>(t)  = 
l^™F(x)dx  which  is  Eq.  (20)  of  Ref.  [5]  for  N  or  M=l.  This  expression  together 

with  the  asymptotic  relation  lim  F(t)  =  VW  lim  (J)(t)  leads  to 

t->oo  t-V® 

0(t)  a  t"®  j  t-*’®  (7b) 

(Eq.  (30)  of  Ref.  [5]),  this  result  is  reiterated  by  Shlesinger  as  Eq.  (21)  in 
a  more  recent  Reference  [lb] .  Thus,  we  see  from  Eqs.  (7a)  and  (7b)  that  two 
difterent  expressions  are  obtained  depending  on  which  of  the  two  computational 
procedures  is  used.  Advocates  of  CTRW  models  then  must  make  clear  which  pro¬ 
cedure  is  actually  appropriate.  Shlesinger  [lb]  attempted  to  show  the  rele¬ 
vance  of  Eq .  (7b)  to  actual  experimental  observations.  He  quoted  a  decay 

which  has  the  form  of  Eq.  (7b)  with  a=0.3  due  to  electron-hole  recombination 
in  sputtered  a-Si:H.  It  should  be  noted  that  Shlesinger' s  calculations  [5] 
leading  to  Eq.  (7b)  as  well  as  the  other  calculations  quoted  in  Sec.  4  of  Ref. 


[lb]  are  based  on  the  assumption  of  unlmolecular  reactions.  On  the  other 
hand,  the  experimentalists  in  Ref.  [25]  of  Ref.  [lb]  have  explicitly  stated 
that  electron-hole  recombination  in  a-Si:H  is  a  bimolecular  process.  This 
conclusion  is  based  on  measurements  of  the  laser  intensity  dependence.  Com¬ 
parison  of  bimolecular  data  with  results  of  a  unimolecular  model  surely  cannot 
be  used  to  resolve  the  dilemma  posed  by  the  two  alternatives  for  computational 
procedure. 

There  is  a  drastic  difference  in  the  long  time  dependences  of  the  relaxa¬ 
tion  function  predicted  by  the  SM  model  depending  on  the  value  of  the  number 
of  defects  per  frozen  dipole  N.  Thus  the  long  time  dependence  of  the  relaxa¬ 
tion  function  obtained  in  the  SM  model  is  not  robust  against  the  choice  of  N. 
Further,  it  does  not  seem  possible  for  SM  to  make  an  estimate  for  the  size  of 
N  in  any  real  glass  or  polymer  system  because  the  identities  of  the  defect  and 
the  frozen  dipoles  have  not  been  specified  for  any  system.  Before  the  SM 
model  can  be  accepted  as  a  general  explanation  for  KWW  relaxation,  validity  of 
the  condition  of  N»1  must  be  checked  for  each  glass  and  polymer  system.  In 
fact,  there  are  many  real  glass  systems  in  which  the  size  of  a  dipole  is  com¬ 
parable  to  that  of  the  smallest  molecular  unit.  In  each  of  these  cases,  it 
would  be  difficult  to  visualize  mobile  defects  which  are  far  more  numerous 
than  dipoles,  and  it  is  not  likely  to  have  N>>1.  The  SM  model  would  then  lead 
to  a  relaxation  function  whose  tail  has  the  dependence  (7)  which  is  not  KWW 
like. 

In  order  to  specify  the  time  scale  for  the  onset  of  the  KWW  behavior,  the 
full  4;(t)  is  needed  and  not  just  its  long-time  behavior.  It  also  should  be 
noted  that  the  SM  model  does  not  derive  the  KWW  function  from  fundamental 
properties  of  dynamical  equations  of  motion.  Within  the  model,  the  power  law 
t®  in  the  KWW  function  Eq .  (6)  is  related  to  an  assumed  power  law  t”^“®  in  the 


pausing  distribution  Eq.  (1).  One  power  law  Is  simply  traded  for  another. 
Even  if  the  N»1  hypothesis  could  be  justified  for  some  systems,  physical 

insight  into  the  KWW  form  would  require  physical  insight  into  the  4;(t) 
function.  Further,  KWW  relaxation  is  observed  in  many  other  relaxation  pro¬ 
cesses  besides  dipolar  relaxations  of  amorphous  polymers  and  glasses.  For 
example,  it  has  been  accurately  established  for  thermoremanent  magnetization 
in  spin  glasses  [11-13]  ,  where  we  can  hardly  think  of  defects  and  diffusion 
processes.  There  seems  to  be  no  basis  for  the  introduction  of  a  pausing  time 
distribution  ();(t)  in  this  case. 

The  discussions  of  this  section  have  revealed  that  the  derived  relaxation 
function  depends  sensitively  on  the  assumption  of  the  form  of  (J;(t)  and  on 
whether  the  condition  of  N»1  is  satisfied.  The  choice  of  the  pausing  time 
distribution  of  the  form  4)(t)~t“^“o  for  long  times  will  be  shown  in  a  later 
section  to  contradict  direct  experimental  evidences  in  the  closely  related 
problem  of  hopping  ions  in  glass  systems.  The  picture  of  mobile  diffusing 
defects  is  very  specific  but  is  immediately  inapplicable  in  the  closely 
related  problem  of  spin  glasses  which  has  been  experimentally  observed  to 
follow  a  full  KWW  form  [11-13]  .  Therefore,  we  conclude  that  the  SM  model  and 
other  CTRW  based  models  are  not  robust  in  their  invocation  of  a  specific  ad 
hoc  pausing  time  distribution,  in  their  requirement  for  a  very  large  ratio  of 
defects  to  relaxors  and  in  the  absense  of  physical  counterparts  for  the  com¬ 
ponents  of  the  models. 

C .  Technical  Results 

In  this  section,  we  shall  follow  the  SM  and  Bendler-Shlesinger  (BS)  [2] 
models  and  work  out  some  of  their  technical  consequences.  From  these  results 
we  will  show  explicitly  that  these  models  can  be  used  to  derive  a  relaxation 
function  which  must  depart  from  the  KWW  form  except  in  the  tail  portion  of  the 
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relaxation  function.  The  early  and  major  part  of  the  derived  relaxation  func¬ 
tion  will  not  have  the  KWW  form  In  disagreement  with  experimental  data. 

In  [1] ,  SM  motivate  the  asymptotic  form  of  the  pausing  time  distribution 
in  Eq.  (1)  by  referring  to  a  similar  hypothesis  used  earlier  [14,15]  in  dis¬ 
persive  transport  in  semiconductors] 16] .  One  of  the  models  used  to  produce  a 
4)(t)  with  a  t"!"®  tail  in  dispersive  tranpsort  was  a  distribution  of  activated 
barriers  [17] .  The  distribution  of  barriers  leads  to  a  distribution  of  rates 
p(\)  which  determines  4)(t)  as 

4^(t)  -  ^  exp(-Xt)p(\)d\ 

Recently,  Bendler  and  Shlesinger  (BS)  [2]  have  advocated  using  this  as  a  model 
within  the  SM  model  to  provide  a  normalizable  pausing  time  distribution  func¬ 
tion.  A  thermally  activated  jump  rate  in  the  presence  of  a  single  free  energy 
barrier  Fq,  XQ=X<aexp(-Fg/kT) ,  is  modified  by  the  entropy  and  energy  fluctua¬ 
tions  generated  by  local  neighboring  rearrangements:  6F  =  6E-T6S.  For  sim¬ 

plicity  the  linear  relation  6S  »  <j6E  is  assumed  by  BS.  As  a  further  assump¬ 
tion  6E  is  taken  to  have  exponential  distribution.  A  generalization  of  the 
distribution  used  by  BS  can  be  written  in  the  form 
exp(-q6E)  ,  0<6E<6E2 

g(6E)  »  A  exp(+q6E)  ,  6Ei<5E<0 

0  ,  6E>6E2,  5E<6Ei  (9) 

Here  A  is  a  normalization  constant.  The  normalization  of  g( 5E)  insures  the 
normalization  of  4)(t)  in  Eq.  (2).  BS  actually  specialized  Eq .  (9)  by  assuming 
6Ei=0  and  6E2*®.  In  other  words,  they  cut  off  the  negative  fluctuations  with¬ 
out  physical  justification.  More  general  distributions  of  the  form  Eq.  (9) 
had  been  studied  over  twenty  years  earlier  by  MacDonald  [18]  in  the  context  of 
modelling  mechanical  relaxations.  Using  p(X)“g( 5E) |d6E/dx|  in  Eq .  (8)  leads 


4;(t)/X  ■  c  [./“^duu”®  exp(-ux)+  r^duu®  exp(-UT)l 

0  a  1  u2 


where 


T=\ot»  a^Tq/(l-cjT), 


ui=exp[-q6Ei/a] , 


U2=exp[-q6E2/ a] , 


and 


C(jj=a/[2-'Ui“a-u2«]  • 

Equation  (10)  provides  an  expression  for  (l;(t)  that  allows  examination  of 
the  properties  of  the  SM  model.  The  BS  choice  6E|^“0  and  6E2*®  corresponds  to 
ui»l  and  U2=“0.  With  this  choice,  the  asymptotic  expansion  of  Eq.  (9)  at  long 
times  yields  4;(T)/X.o'^“r(l+a)i:"^“®,  which  is  a  power  law  tail  as  required  by 
the  SM  model  (as  long  as  a*kTq/(l-aT)<l) .  As  discussed  in  Section  B,  this 

leads  only  to  the  tall  end  of  the  KWW  function  of  Eq.  (6)  with 

d  -  [0.659c/ar(a)r(l-o)r(l+o)]\“  (11) 

To  examine  the  timescale  for  the  onset  of  the  tail  of  the  4i(t)  given 

by  Eq.  (10)  we  have  evaluated  Eq .  (10)  numerically  for  several  values  of  a. 
Fig.  1  shows  the  calculated  4)(t)/\o  for  a=0.5  with  the  BS  choice  u^^l,  U2“0. 
The  dashed  curve  shown  is  the  asymptotic  expression  ar(l+a) which 


approaches  the  calculated  4)(i:)/\o  values  (solid  curves)  near  As  a  meas¬ 

ure  of  the  meaning  of  this  onset  time,  we  have  also  computed  the  median  time 


Tai,  defined  by 

y  dT(|^(T)/X 

CO  o 

The  median  time  divides  the  total  weight  of  the  (^(t) 


(12) 

distribution,  giving 


exactly  half  the  weight  for  the  short  time  region  ^<^31  and  half  the  weight  for 


the  long  time  region  T>Tn.  The  curve  in  Fig.  2  is  shaded  for  times  up  to  Tm- 
Notice  that  (1;(t:)  assumes  the  form  only  beyond  ^^1.  Since  the  SM  model 


by  SM.  The  experimentally  measured  quantity  Is  the  frequency  dependent  die¬ 
lectric  function  e*((i))*e' (u)~le"(u) .  In  the  theory  of  dielectric  relaxation: 

e*(u))-^rdt  e"^“‘'[-d4)/dt] 

The  polyvinylacetate  experimental  e*(u))  data  Is  well  fit  by  Eq.  (13)  with  a 
KWW  function  <t)  of  Eq.  (6)  having  the  (slightly  temperature  dependent)  frac¬ 
tional  exponent  a~0.56  [19].  Moreover,  the  fit  covers  significant  frequency 
ranges  on  both  sides  of  the  frequency  ojp  of  the  dielectric  loss  peak.  The 
experimental  frequency  range  Is  wide  enough  to  cover  essentially  the  entire 
relaxation  process  from  beginning  to  end.  However,  the  SM  model  will  not 
produce  a  KWW  function  for  T<Tni  and  therefore  the  model  will  not  fit  the 
dielectric  £'((jj),  e'*(oj)  data  at  least  for  the  high  frequency  side  of  cop.  This 
can  be  understood  from  the  fact  that  ojp  for  the  KWW  function  corresponds 
approximately  to  the  time  at  which  the  decay  has  reached  1/e  of  Its  Initial 
value.  This  disagreement  of  the  SM  result  for  the  high  frequency  side  of  the 
dielectric  relaxation  function  also  occurs  for  the  primary  relaxations  of 
supercooled  liquids  and  polymers  stiidled  by  Williams  and  Watts  themselves 
[9,20] . 

Dielectric  relaxation  can  be  carried  out  In  the  time  domain  by  measurement 
of  discharge  of  the  depolarization  current  for  a  wide  range  of  dielectric 
materials.  The  current  In  the  short  time  domain  follows  the  power  law  of  time 
dependence. 

i(t)  a  t"n  (14) 

This  form  of  the  dielectric  response  has  been  known  since  the  turn  of  the 
nineteenth  century  under  the  name  of  the  Curie-von  Schweidler  law  [21,22].  At 
longer  times,  the  time  dependence  of  the  current  departs  from  Eq.  (14).  From 
the  formalism  of  dielectric  relaxation,  i(t)  a  d())/dt.  The  classic  Curie-von 
Schweidler  law  is  consistent  with  the  KWW  form  Eq.  (6),  because  d(})/dt* 


I 


I 
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“acit®''^exp(-cit®)  and  at  short  times  the  t®"^  factor  dominates.  This  demon¬ 
strates  how  essential  It  Is  for  any  dielectric  relaxation  model  to  capture  the 
entire  KWW  form  over  both  the  short  and  long  time  regions.  Zwanzlg  [23]  has 
emphasized  this  In  his  comment  on  another  model. 

It  Is  Interesting  also  to  consider  the  effect  of  modifying  the  choice  for 
g(6E)  In  the  BS  model  on  the  median  time.  We  have  calculated  Eq.  (10)  numeri¬ 
cally  for  the  distributions  Eq.  (9)  with  6E2*®  and  various  values  of  (5Ei<0, 
which  corresponds  to  U2“0  and  ui>l.  Including  contributions  with  6Ei<0  fur¬ 
ther  delays  the  onset  of  the  tall  beyond  the  median  time.  Fig.  2  shows 
the  calculated  4;('c)/X<j  corresponding  to  the  case  where  the  fraction  of  the 
weight  of  g(6E)  with  6E<0  Is  f_-0.45  (l.e.  nearly  symmetric  In  positive  and 
negative  fluctuations).  The  asymptotic  power  law  tall  Is  again  shewn  by  the 
dashed  curve.  The  onset  of  c^Ct)  to  the  power  law  tall  Is  delayed  even  further 
beyond  the  median  time  xn  as  compared  to  the  f__=»0.0  case  In  Fig.  1.  As  shown 
in  Fig.  3,  the  value  of  tm  in  fact  decreases  rapidly  as  the  6E<0  contributions 
occupy  more  weight  in  the  g(5E)  distribution,  l.e.  as  f__  increases.  However, 
they  do  not  contribute  to  the  tail  except  in  the  normalization  prefac¬ 
tor;  (|;(x)-*-ar(l+a)(2-f_)x“^“®.  This  modified  prefactor  moves  the  onset  of  the 
power  law  tail  to  shorter  times,  but  not  as  rapidly  as  the  decrease  ot  x^. 
The  result  is  that  as  the  fluctuations  5E  in  g(6E)  are  allowed  to  becocffi  more 
symmetrical,  the  SM  model  produces  a  relaxation  function  that  coincides  with  a 
KWW  form  over  a  steadily  diminishing  tail  portion  of  the  dielectric  relaxation 
curve.  This  situation  is  further  exacerbated  if  the  positive  fluctuations  are 
truncated  at  a  finite  value  of  6E2,  and  also  if  a  Is  Increased.  Thus  the  SM 
model  with  (i;(t)  in  turn  modelled  by  BS,  even  aside  from  the  status  of  its 
assumptions,  can  only  lead  to  a  relaxation  function  that  coincides  with  the 
KWW  form  in  a  long  time  regime  that  cannot  coincide  with  the  full  time  regime 


in  which  the  observable  relaxation  occurs. 

13 


As  seen  in  Eqs.  (6)  and  (11),  the  KWW  relaxation  time  t*,  defined  by  re¬ 
writing  the  KWW  form  as  4i(t)-exp[-(t/i:*)“]  is  calculated  by  BS  to  be 
BS  state  that  (J<(t)/Xo  goes  to  a  constant,  (l;(o) /Xo“a/(l+a) ,  for  t<1*  If 


(|;(t)=a/\o/(l+cx)  at  short  times,  then  Eq.  (5)  leads  to  I(u) (!■*■«) 
large  u  for  any  lattice.  This  prochices  a  short  time  exponential  decay  as 
claimed  by  BS.  In  particular,  (|)(t)=exp(-t/To)  with  To“(l+a)Xo“^/ac*  This 
would  lead  to  the  prediction  t*  a  Tq,  which  disagrees  with  a  wide  range  of 
experimental  data  [24-26]  .  Note  from  Fig.  1  that  (j<('c)/Xo  is  really  not  con¬ 
stant  for  T<1  as  stated  by  BS,  but  Instead  decays  by  almost  60%  of  its  initial 
value  during  this  interval.  The  initial  decay  of  (Ji(T)/\o  becomes  even  steeper 
if  5E<0  is  allowed  as  in  Fig.  2.  Therefore,  the  short  time  behavior  predicted 
by  BS  model  will  not  actually  be  an  exponential  decay  for  the  g(6E)  distribu¬ 
tion  as  chosen  by  them. 

The  unsatisfactory  nature  of  the  SM  model  with  4<(t)  as  modelled  by  BS  in 
deriving  only  a  tail  portion  of  the  KWW  form  is  not  unique  to  the  activated 
barrier  model  for  (|<(t).  Any  (l;(t)  with  an  asymptotic  tall  will  produce 

the  same  problems  as  discussed  above.  For  example,  in  an  earlier  discussion 
in  the  context  of  dispersive  transport,  Shleslnger  [15]  adopted  a  Levy  stable 
function  in  time  (see  Appendix  A  of  Ref.  [15])  for  4;(t).  This  has  the  pro¬ 
perty  that  the  Laplace  transform  of  4;(t)  is  a  fractional  exponential  in  the 
Laplace  transform  variable  u;  L[4,(t)  ]»exp[-(u/\o)a] .  Here  Xq  is  a  rate  para¬ 
meter  that  is  distinct  from  that  used  in  the  barrier  distribution  model.  For 
example,  when  a“0.5,  the  corresponding  form  of  (J)(t)  can  be  written  analyti¬ 
cally:  (i)(t)  =•  exp[-l/4Xot] /  2>^x( Xot)'^.  This  is  properly  normalized  and  has  a 
long-time  t"^"®  tail.  Indeed,  the  only  physical  rationale  implied  for  the 
choice  of  the  Levy  stable  function  is  that  it  has  the  appropriate  long-time 
tall.  We  have  calculated  this  4)(t)  with  a“0.5  and  the  results  are  shown  in 
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Fig.  4.  Again,  the  onset  of  the  tail  is  found  to  be  beyond  the  median 


time  Zfn-  Note  that  a  Levy  stable  function  in  time  for  ci;(t)  is  quite  different 


from  the  use  of  stability  for  the  frequency  spectrum  of  (fC^)  discussed  by 


us  elsewhere  [27]  .  Recall  that  we  have  shown  that  the  relaxation  functions 


derived  by  SM  and  BS  cannot  e3q>laln  the  short  time  regime  of  the  KWW  function 


which  is  observed  so  widely  that  it  is  called  the  Curie-von  Schweidl^^sr  law. 


For  the  same  reason  one  might  question  also  the  claims  in  Refs.  [14,15]  that 


the  characteristic  dependence  of  the  transient  current  in  dispersive 


transient  transport  in  chalcogenlde  glasses  [17]  for  t<<tx,  the  transit  time. 


can  be  explained  by  a  (j)(t)~t“^“®  at  long  times  in  CTRW  framework. 


In  this  section,  the  technical  results  were  carefully  carried  out  and 


examined  in  detail.  We  found  that  the  result  of  all  the  modelling  of  (j;(t) 


[1,2]  will  produce  relaxation  functions  which  coincides  with  the  KWW  only  over 


a  tall  portion  at  long  times.  This  contradicts  the  dielectric  relaxation  data 


in  the  time  domain  (Curie-von  Schweidler  law)  as  well  as  the  frequency  depen¬ 


dence  of  the  measured  dielectric  function  c*(a))  first  given  by  Williams  and 


Watts  themselves  [9,20] .  We  conclude  that  the  technical  results  of  the  SM  and 


the  BS  or  any  of  the  presently  proposed  models  for  ii(t)  do  not  produce  the 


entire  KWW  form.  Thus,  CTRW  models  of  relaxation  do  not  provide  a  satisfac¬ 


tory  basis  for  description  of  physical  phenomena. 


D.  Physical  Tests  of  Assumptions 


One  of  the  key  assumptions  of  the  SM  model  is  that  the  pausing  time  dis- 


trilxition  has  a  t“^”®  tail.  The  repeated  advocation  of  the  t”^“®  form  [1,2, 


14,15]  gives  the  impression  that  it  is  generally  applicable  to  different  ran¬ 


dom  walkers  in  many  different  systems.  We  have  found  physical  systems  which 


are  glass-forming  materials  and  have  mobile  random  walkers.  The  pauslng-tlme 
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dlstributloa  4»(t)  applies  directly  [28,29]  in  these  systems,  if  the  experimen¬ 
tal  data  is  Interpreted  in  a  CTRW  framework  [3] .  However,  we  must  emphasize 
that  a  CTRW  framework  is  not  required  to  Interpret  the  experimental  data. 
Nonetheless,  if  a  CTRW  framework  is  used,  it  can  be  put  to  a  clean  test  in 
such  systems. 

To  do  this  we  examine  electrical  relaxations  due  to  ionic  diffusion  in 
glasses  above  and  below  the  glass  transition  temperature.  It  was  recognized 
by  Moynlhan,  Macedo  and  co-workers  [28,29]  that  for  electrical  relaxations  in 
dielectrics  containing  a  substantial  concentration  of  mobile  charges,  the 
appropriate  physical  picture  is  the  decay  of  the  electric  field  ^  at  constant 
displacement  vector  In  contrast,  relaxation  of  permanent  dipoles  is  char¬ 
acterized  by  the  decay  of  D  under  the  condition  of  constant  The  decay  of 

the  electric  field  is  due  to  diffusion  or  migration  of  the  mobile  ions  and  may 
be  described  by 

E(t)  -  E(o)Y(t)  (15) 

where  y(t)  is  a  decay  function.  If  the  mobile  ions  are  described  as  random 
walkers  in  the  continuous  time  random  walk  formalism,  then  ({;(t)>-dY(t) /dt. 
That  is,  the  rate  of  decay  of  the  electric  field  is  a  measure  of  the  pausing 
time  distribution  of  the  mobile  ions.  Experimentally,  systems  including 
alkali  network  oxide  glasses,  fused  salts,  and  concentrated  aqueous  electro¬ 
lyte  solutions  have  been  studied  for  electric  field  relaxation  [27-32] . 

Let  us  focus  our  attention  on  the  alkali  silicate  glasses.  There  are  no 
permanent  dipoles  in  these  systems  and  dielectric  relaxation  is  due  to  diffu¬ 
sion  of  the  alkali  ions.  Since  there  are  no  frozen  dipoles  in  an  alkali  sili¬ 
cate  glass  and  only  a  diffusing  species,  we  have  a  direct  test  of  the 
(|;(t)  a  t”^“®  hypothesis  of  the  SM  model.  We  need  not  worry  about  the  validity 
of  the  subsidiary  hypotheses  of  the  model.  If  we  take  the  SM  model  seriously. 


in  particular  that  the  t"l”a  pausing-time  distribution  is  very  general,  then 
the  diffusion  of  the  alkali  ions  should  be  well  described  by  Eq.  (1).  How¬ 
ever,  experimentally  it  has  been  found  repeatedly  [28-32]  that  the  decay 
function  -dY(t)/dt  in  Eq.  (14)  is  not  t“l“o  but  remarkably  is  itself  the  time 
derivative  of  a  KWW  fractional  exponential  function.  Examples  include  alkali 
silicate  and  borate  glasses,  a  glass-forming  40mol%  Ca(H03)2-60mol%  KNO3  melt, 
Na  beta  alumina,  LlCl  electrolyte  solutions  and  others  [28-32,24].  These  are 
examples  of  KWW  decay  involving  only  diffusing  species.  This  means  that  if  we 
wish  to  use  a  pausing  time  distribution  to  describe  these  experiments,  it  must 
have  the  form  (j)(t)“-dY(t)/dt*( atot“J-/'t*®)exp(-[ t/x*]®) •  This  obviously  does 
not  have  a  t~^““  long  time  tall.  Thus  the  claim  of  generality  of  a  t“^”®  tail 
in  the  pauslng-tlme  distribution  is  not  supported  by  experimental  data  on 
electric  field  relaxations  due  to  diffusing  ions  in  glasses. 

The  failure  of  a  (l)(t)  with  a  tail  to  describe  experimentally  the 

diffusing  ions  can  also  be  seen  explicitly  by  calculation  of  a  corresponding 
frequency  spectrum  known  as  the  electric  modulus,  M*(u).  This  is  the  inverse 
of  the  dielectric  function,  M*(a))=l/e*(u)),  and  is  related  to  4;(t)  by 

M*(oj)  -  M’(a))+iM"((j)  -  M,[l-^/»dte"^"Vt)]  (16) 

Experimentally  [28-32]  ,  M"((jj)  shows  a  characteristic  asymmetric  broad  peak 
which  is  skewed  toward  high  frequencies.  The  data  can  be  fit  very  well  if 
(l>(t)  is  taken  to  be  the  derivative  of  a  KWW  function  [28-32]  .  Now  let  us  com¬ 
pare  this  with  the  results  of  the  SM  model  with  (l)(t)  as  modelled  by  BS.  Using 
the  activated  barrier  distribution  model  for  (j;(t),  Eq .  (10),  we  find 

M'(u))/M„  -  4  d?) 

M"(aj)/M„  -  c^Q  [i/'^du  +  ^^^du  -H_] 

where  QHu/Xq.  We  have  evaluated  Eqs.  (17)  and  (18)  numerically  for  different 


choices  of  the  distribution  Eq .  (9).  The  calculated  M"(oj)  shows  the  peak  near 


oj^Xq  shape  of  the  curve  disagrees  with  experiment.  If  the  fluctua¬ 
tions  6E  are  symmetric  In  positive  and  negative  contributions  (l.e. ,  ui»u2“^) 
the  calculated  M"((i))  curve  Is  symmetric  about  Its  peak.  This  Is  shown  In 
Fig.  5.  On  the  other  hand  the  experimental  M"  curves  are  skewed  as  shown  In 
Fig.  7.  If  we  follow  BS  and  remove  all  negative  fluctuations,  emphasizing  the 
power  law  tall  as  much  as  possible,  then  the  calculated  M"(u3)  curve  Is  asym¬ 
metric  as  shown  In  Fig.  6  but  skewed  In  the  opposite  way  from  the  experimental 
curves.  Hence  the  BS  model  of  4;(t)  based  on  a  distribution  of  barriers  with 
only  positive  fluctuations  as  well  as  more  physically  reasonable  models  of 
(j;(t)  which  have  symmetric  fluctuations  all  disagree  with  the  experimental 
measurements.  Therefore,  diffusing  species  In  glass  forming  systems  are  not 
generally  described  by  a  t“^“®  tall.  The  hypothesis  of  the  SM  model  that  (i;(t) 
of  a  yet  unspecified  defect  has  a  t“^“«  tall  thus  becomes  questionable. 

There  Is  another  difficulty  in  using  a  defect  diffusion  model  for  relaxa¬ 
tion  In  realistic  glass  forming  materials.  In  many  glasses  and  amorphous 
polymers  [33] ,  several  relaxations  of  different  nature  are  present  In  the  same 
system.  It  has  been  demonstrated  that  even  small  molecule  glasses  have  both 
primary  and  secondary  relaxations  [34]  .  Often  the  different  relaxations  in 
the  same  system  all  have  the  KWW  form,  although  each  has  its  own  own  charac¬ 
teristic  fractional  exponent  a  [35]  .  This  is  known  to  proponents  of  the  CTRW 
based  relaxation  models.  In  fact  they  [36]  have  fit  the  experimental  data  of 
a  secondary  relaxation  of  a  small  molecule  glass  to  the  KWW  form  and  claimed 
good  agreement.  Consider,  for  example,  glassy  polycarbonate.  There  are  at 
least  three  distinct  relaxations  traditionally  labelled  alpha,  beta,  and  gamma 
[37,33].  The  alpha  and  gamma  relaxations  have  been  extensively  studied 
[33,  35,37,38],  and  they  both  have  the  KWW  form  but  with  quite  different 
values  for  the  fractional  exponent:  a=0.4  for  the  alpha  relaxation  [35]  and 


a=K).2  for  Che  gamma  relaxation  [38] •  Many  other  such  examples  can  be  cited. 
This  type  of  example  presents  a  difficulty  for  the  SM  model,  because  the  value 
of  a  in  the  VUV  tail  derived  is  the  same  as  that  in  the  hypothesized  t"l“® 
tail  of  the  pausing-time  distribution  of  the  diffusing  defect,  Eq.  (1).  In 
order  to  account  for  two  different  KWW  exponents  observed  in  the  same  physical 
system,  the  SM  model  would  require  two  different  types  of  diffusing  defect, 
each  with  a  different  t~l~®  tail.  Furthermore,  defect  1  (2)  must  be  required 
to  interact  only  with  dipoles  of  the  alpha  (gamma)  relaxation  but  never  with 
the  gamma  (alpha)  relaxation.  Such  a  requirement  seems  very  difficult  to  jus¬ 
tify  on  physical  grounds. 

E.  Predictability 

An  important  aspect  of  predictability  is  whether  significant  phenomena 
that  occur  empirically  are  predicted  or  at  least  included  in  the  phenomenology 
of  Che  model.  We  have  seen  the  models  of  relaxation  based  on  CTRW  have  not 
predicted  an  entire  KWW  function  but  rather  a  relaxation  function  that  ap¬ 
proaches  a  KWW  form  exp(-cit®)  only  at  the  tail  end.  Therefore,  they  [1,2] 
cannot  even  claim  that  their  models  predict  a  genuine  KWW  relaxation  function. 
One  of  Che  assumptions  Chat  has  to  be  made  is  Che  pausing  time  distribution 
4;(t)  having  a  long  time  Call  of  t"^“®.  Since  the  KWW  relaxation  function  is 
generally  observed,  a  fact  subscribed  by  authors  of  these  CTRW  models,  logi¬ 
cally  the  key  assvanption  of  4;(t)~t“l”i  must  be  at  least  or  even  more  ubiqui¬ 
tous.  The  electrical  relaxation  data  of  many  glass-forming  systems  (Section 
D)  demonstrate  that  the  latter  is  not  true.  Hence  this  key  assumption  is 


questionable.  Thus,  at  the  expense  of  at  least  one  ad  hoc  assumption,  the  SM 
model  has  predicted  only  part  of  a  KWW  relaxation  function.  There  are  no 
other  concrete  predictions  or  further  applications. 


The  BS  model,  with  the  assumption  of  a  one-sided  fluctuation  of  energy 
barrier  heights  with  an  exponential  distribution  does  provide  an  additional 
prediction.  There,  the  quantity  a  In  the  t”l"“  tall  of  4)(t)  becomes  tempera¬ 
ture  dependent  and  the  exponent  a  of  the  predicted  tall  of  the  KWW  function, 
exp(-cit®)  will  also  be  temperature  dependent.  Let  us  compare  this  prediction 
with  experimental  data.  Actaully  there  Is  no  completely  consistent  way  of 
making  such  a  comparison.  The  experimental  data  show  that  the  relaxations  are 
entire  KWW  functions,  while  the  BS  model  predicts  only  relaxation  functions 
whose  tails  are  KWW  like  (see  Section  C).  Even  when  we  Ignore  this  problem  of 
the  BS  model,  experimentally  there  are  glass-forming  materials  such  as  SI-O2 
[39],  polystyrene  [40],  polyethyl  acrylate  [41],  etc.  for  which  the  fractional 
exponent  a  of  the  entire  KWW  relaxation  function  In  each  case  remains  constant 
over  a  sizeable  experimental  temperature  range.  The  polyethyl  acrylate  data 
Is  particularly  worth  noting  because  It  was  obtained  by  Williams  and  Watts 
themselves.  McDuffie,  Quinn  and  Litovltz  [42]  have  studied  the  dielectric 
properties  of  glycerol-water  mixtures.  They  observed  that  the  shape  of  the 
relaxation  function  Is  Independent  of  the  concentration  and  conclude  this 
result  Is  not  In  accord  with  the  concept  that  the  origin  of  the  departure  from 
exponential  relaxation  lies  In  the  existence  of  a  distribution  of  barrier 
heights  in  the  liquid.  Litovltz  and  >fcDuffie  [42]  compare  the  dielectric 
rotational  time  with  structural  break-up  time  and  point  out  a  difficulty  of 
the  Glarum's  model  [7]  when  applied  to  associated  liquids.  There  are  also 
many  other  systems  whose  relaxation  behavior  are  of  the  KWW  form,  with  tem¬ 
perature  dependent  fractional  exponents  a  [43]  .  However,  again  in  each  of 
these  cases,  the  relaxation  can  be  described  by  an  entire  KWW  function  In  dis¬ 
agreement  with  the  predicted  relaxation  function  of  the  BS  model. 


There  are  three  empirical  relations  of  relaxation  in  addition  to  the  KWW 


function  (which  we  call  the  first  relation)  that  have  been  seen  repeatedly  in 
experimental  data.  The  second  relates  the  effective  relaxation  time  x*  that 
appears  in  the  KWW  form,  now  written  as 

0(t)  -  exp(-(t/x*)a)  ,  (19) 

to  a  primitive  relaxation  time  xq  that  has  a  microscopic  and  fundamental 
meaning  as 

.  1-a  .1/a 

“  ^“c  ■'o^  (20) 

where  uc  is  a  characteristic  frequency.  It  is  important  to  emphasize  that  the 
fractional  exponent  a  and  x*  appearing  in  Eq.  (20)  are  the  same  as  those  in 
the  KWW  form  (19).  The  third  relation  states  that  for  the  same  relaxation 
process,  the  relaxation  function  is  an  entire  KWW  function  (19)  if  i^xo>>l  but 
crosses  over  to  an  entire  linear  exponential 

<t)(t)  »  exp(-t/xo)  (21) 

if  a)cTo<<l.  The  quantity  xq  in  Eq.  (21)  and  cog  in  the  conditions 

again  the  same  as  that  appears  in  Eq.  (20).  The  relaxation  models 
based  on  CTRW  predict  neither  of  these  two  extra  empirically  verified  rela¬ 
tions  [24,26,30,31,35,38,44,45]  but  at  best,  predicts  x*  a  Tq  [25]. 


F .  Summary  and  Ck^nclusions 

We  have  initiated  an  examination  of  existing  models  of  relaxation  by  con¬ 
sideration  of  models  [1,2]  based  on  the  continuous  time  random  walk  (CTRW) 
formalism  [3] .  CTRW  is  only  a  mathematical  scheme.  To  address  a  physical 
problem  such  as  relaxation  in  glass  forming  systems  described  by  the 
Kohlrausch-Williams  and  Watts  or  KWW  function,  additional  physical  inputs  are 
required.  These  include  (1)  the  pausing  time  distribution  cli(t),  (2)  the 
ratio  of  defects  to  relaxors,  and  (3)  the  contact  nature  of  the  def ect-relaxor 
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interaction  that  causes  relaxation.  In  order  to  obtain  a  physical  model  of 
relaxation  based  on  CTRW,  these  inputs  must  be  the  consequences  of  or  deriv¬ 
able  from  some  general  physical  principles  and  fundamental  laws.  This  level 
is  not  achieved  in  either  of  the  two  proposed  models  [1,2].  In  the  proto¬ 
typical  work  by  Shlesinger  and  Montroll  (SM)  [1] ,  the  choice  made  for  the 
pausing  time  distribution  function  (l;(t)  has  a  t“l~o  tail  (o<a<l)  at  long  times 
is  made.  The  exact  nature  of  (lj(t)  at  earlier  times  before  the  t“l~<i  tail 

takes  over  is  left  unspecified.  An  unspecified  contact  interaction  between 

defect  and  relaxor  is  assumed.  Another  essential  assumption  is  that  the  num¬ 
ber  of  defects  per  relaxor  is  large.  With  these  hypotheses,  SM  have  shown  us 
how  to  turn  the  crank  of  CTRW  mathematics  to  derive  a  tail  portion  of  the  KWW 
relaxation  function.  This  result  depends  on  how  I(u)  in  Eq .  (5)  is  evaluated 
at  small  u.  If  an  earlier  method  for  this  evaluation  due  to  Shlesinger  [5]  is 
used,  the  result  would  be  given  by  relation  (6S),  i.e.  4)(t)~exp(-C2t''“)  for 
and  it  would  not  be  KWW  anywhere.  The  CTRW  workers  [1,5]  should  clarify 
the  methodology  to  be  employed  in  evaluating  I(u)  for  small  u. 

Even  with  the  machinery  used  by  SM,  one  only  obtains  a  relaxation  function 
that  coincides  with  the  KWW  function  in  a  limited  long  time  regime.  The 
shorter  time  regime  in  which  the  KWW  form  is  empirically  satisfied  is  not 
described  adequately.  Thus  Eq .  (29)  of  the  SM  paper  is  imprecise.  It  is  an 

asymptotic  equation  for  the  behavior  of  4>(t)  for  large  t  and  should  not  be 

interpreted  as  an  equation  valid  for  all  times.  To  make  sure  that  this  dis¬ 
tinction  between  an  asymptotic  relation  and  an  equation  is  not  splitting 
hairs,  we  examined  the  two  pausing  time  distributions  respectively  modelled  by 
Bendler  and  Shlesinger  (BS)  [2]  and  by  Shlesinger  above,  [15] .  Both  of  these 
are  distributions  which  have  a  t”l~a  tail  and  normalized  over  all  times.  We 
found  that  the  weight  of  the  distribution  for  the  short  time  regime  where 
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departs  strongly  from  the  t“^“®  dependence  has  more  than  half  of  the  total 
weight  of  the  distribution.  These  calculations  point  to  a  serious  problem 
common  to  all  existing  models  of  relaxation  based  on  CTRW.  They  have  not 
accomplished  the  primary  objective  of  a  derivation  of  an  entire  KWW  relaxation 
function.  Their  relaxation  function  will  approach  a  KWW  form  only  at  times 
after  the  quantity  to  be  relaxed  has  been  reduced  by  more  than  one-half  or 
much  more  than  half  of  Its  Initial  unrelaxed  value.  Experimentally  this 
Initial  relaxation  as  well  as  Its  continuation  to  long  times  correspond 
throughout  to  a  single  KWW  function.  This  Initial  part  of  the  KWW  relaxation 
Is  observed  repeatedly  In  dielectric  relaxation  time  domain  measurements  based 
on  the  technique  of  charging  and  discharging  currents.  The  phenomenon  has 
been  referred  to  as  the  Curle-von  Schweldler  [21]  law.  In  the  frequency 
domain,  the  same  phenomenon  la  observed  by  the  wide  occurrence  of  the 
frequency  dependence  (o<n<l)  for  the  real  and  imaginary  parts  of  the  dielec¬ 
tric  function  e*(a))  on  the  high  frequency  side  of  the  dielectric  loss  peak 
frequency  up.  This  feature  has  been  emphasized  by  Jonscher  [22,46].  It  was 
pointed  out  [35]  that  this  feature  not  only  follows  from  the  KWW  function  but 
also  from  other  empirical  dielectric  relaxation  functions  such  as  those  pro¬ 
posed  by  Cole-Cole,  Cole-Davidson  and  by  Havrlliak-Nagami[45] .  The  relaxation 
functions  predicted  by  models  of  relaxation  [1,2]  based  on  CTRW  contradicts 
this  well  established  phenomenology  of  dielectric  relaxation. 

After  the  assumption  of  (j;(t)~t"l“o  has  been  made,  the  relaxation  function 
obtained  by  CTRW  depends  sensitively  on  whether  number  of  defects  N  per 
relaxor  is  large.  We  have  shown  explicitly  how  the  0(t)  deviates  from  the  KWW 
form  even  In  a  long  time  regime  when  N  Is  finite.  (See  discussions  following 

(7)).  Further,  we  have  pointed  out  possible  physical  cases  where  the 
assumption  of  very  large  N  Is  untenable.  This  lack  of  robustness  of  CTRW 
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models  against  choice  of  N  Is  not  appropriate  for  models  designed  to  describe 
a  very  general  phenomenon,  namely  KUW  relaxation. 

We  have  also  given  evidence  to  dispel  the  belief  that  the  distribution 
with  a  tall  have  general  applicability.  The  evidence  comes  directly 
from  experimental  measurements  of  electrical  relaxation  of  diffusing  Ions  In 
glass  forming  systems.  The  experimental  data  [28-32]  clearly  Indicate  that, 
If  we  Insist  on  a  CTRW  approach,  the  pausing  time  distribution  function  Is  the 
derivative  with  respect  to  time  of  a  KWW  function,  and  this  does  not  have  a 
fl-a  tall. 

Furthermore,  the  CTRW  models  [1,2]  lead  to  a  physically  untenable  picture 
when  applied  to  realistic  glass  forming  materials  which  have  primary  and  sec¬ 
ondary  relaxations  and  both  are  well  described  by  KWW  functions.  The  differ¬ 
ent  fractional  exponents  now  require  two  types  of  defects  each  with  Its  own 
inverse  power  tail.  One  type  of  defect  is  responsible  for  only  one  type  of 
relaxation  process  and  not  the  other.  The  number  of  assumptions  and  hypothe¬ 
ses  involved  increases  rapidly,  making  the  applicability  of  the  models  [1,2] 
less  likely. 

All  theoretical  models  in  physics  are  ultimately  judged  by  the  capability 
to  generate  new  testable  predictions.  The  models  [1,2]  based  on  CTRW  at  the 
expense  of  an  ad  hoc  assumption  of  the  pausing  time  distribution  function  pre¬ 
dicts  a  relaxation  function  that  only  approaches  a  KWW  function  at  long  times. 
At  an  initial  but  experimentally  significant  part  of  the  relaxation,  this  pre¬ 
dicted  relaxation  function  falls  repeatedly  when  compared  with  dielectric 
relaxation  phenomena.  There  Is  no  other  sustainable  prediction.  In  the  model 
by  Bendler  and  Shlesinger  [2]  for  (i;(t)  there  Is  one  additional  prediction. 
The  extra  prediction  is  that  the  fractional  exponent  a  will  necessarily  be  a 
function  of  temperature.  However,  there  are  a  number  of  glass  forming  systems 


in  which  the  experimentally  measured  a  is  insensitive  to  temperature,  in  dis¬ 
agreement  with  this  prediction. 

In  brief,  we  have  examined  the  models  of  relaxation  [1,2]  based  on  CTRW. 
The  conclusion  we  reach  is  that  they  taist  be  rejected  as  relevant  models  for 
KWW  relaxations  in  glass  forming  systems  because  they  lack  robustness,  employ 
certain  questionable  technical  procedures,  introduce  general  assumptions  that 
fail  to  meet  specific  physical  tests,  lead  to  unverified  predictions,  and  fail 
to  predict  or  describe  verified  properties  of  KWW  relaxation. 
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Fig.  1  Solid  curve  is  calculated  from  Eq.  (10)  with  a*0.5  corresponding  to 
g(6E)  in  Eq.  (9)  with  6E2“«  and  only  positive  fluctuations:  f_*0.0. 

Dashed  curve  is  Che  asymptotic  power  law  tall,  xm  is  the  median  time, 
Eq.  (12)  .  Here  x*\oC* 


h*  0.4 


0.2 
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Fig.  2  Solid  curve  is  calculated  from  Eq.  (10)  with  a“0.5,  corresponding  to 
g(5E)  in  Eq.  (9)  with  6E2“«>  a  fraction  of  negative  fluctuations, 
f_*0.45.  Dashed  curve  is  the  asymptotic  power  law  tail.  xm  is  the 
median  time,  Eq.  (12)  .  Here 
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Median  time,  Eq.  (12),  calculated  from  Eq.(lO)  with  a“0.5  correspon¬ 
ding  to  g(6E)  in  Eq .  (9)  with  6E2“®  and  a  fraction  of  negative  fluc¬ 
tuations  f  . 


Fig. 


log(aj/Xo) 

Beal  and  Imaginary  parts  of  the  electric  modulus,  calculated  from  Eqs 
(17)  and  (18)  with  a*0.5  corresponding  to  g(6E)  in  Eq.  (9)  with  sym¬ 
metric  positive  and  negative  fluctuations  6E  (ui»u2“^*2.0) . 
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Iog(w/Xo) 

Fig.  6  Real  and  imaginary  parts  of  the  electric  modulus,  calculated  from  Eq 
(17)  and  (18)  with  a“0.5  corresponding  to  g(6E)  in  Eq.  (9)  with  6E2 
and  only  positive  fluctuations:  f  *0.0. 
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log  f  CHz:> 

M"  vs  log  f  data  at  several  temperatures  for  K20-3Si02  glass  annealed 
at  479“ C  and  quenched.  The  points  are  experimental  data  and  curves 
are  recalculated  by  us  using  Eq.  (16)  and  (i;(t)-exp[-(t/T)«]  with  the 
parameters  a.  Mg,  and  x  given  by  Boesch  and  Moynlhan  (L.P.  Boesch  and 
C.T.  Moynlhan,  J.  Non.  Cryst.  Solids  17 ,  44  (1975)  ). 
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